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Wo study B° ^ p°p° decays in a sample of 465 x 10*' r(4S) BB events eoUeeted with the 
BABAR detector at the PEP-II asymmetric-energy e+e^ collider located at the Stanford Linear 
Accelerator Center (SLAG). We measure the branching fraction B = (0.92 ± 0.32 ± 0.14) x 10"" and 
longitudinal polaxization fraction /i, = 0.75lo^j4 ±0.04, where the first uncertainty is statistical and 
the second is systematic. The evidence for the B° p"p" signal has a significance of 3.1 standard 
deviations, including systematic uncertainties. We investigate the proper-time dependence of the 
longitudinal component in the decay and measure the CP-violating coeflicients Sl" = (0.3±0.7±0.2) 
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and C£° = (0.2 ± 0.8 ± 0.3). We study the implication of these results for the unitarity triangle 
angle a. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



Measurements of CP-violating asymmetries test the 
flavor structure of the standard model by over- 
constraining the Cabibbo-Kobayashi-Maskawa (CKM) 
quark- mixing matrix [l[ . The time-dependent CP asym- 
metry in the decays of or mesons to a CP eigen- 
state dominated by the tree-level amplitude b — > uud 
measures sin2acff, where a^s differs from the CKM uni- 
tarity triangle angle a = arg [-V^^V^^^/V^^V*^] by a quan- 
tity Aa accounting for the contributions from loop (pen- 
guin) amplitudes. The value of Aa can be extracted from 
an analysis of the full set of isospin-related channels . 

Since the tree contribution to the B^ —>■ p^p^ 0] de- 
cay is color-suppressed, the decay rate is much smaller 
than those of the 5° — > p^p^ and chan- 
nels 0, 0, H, 0, Sl> ^'^d is sensitive to the penguin 
amplitude. Therefore a stringent limit on Aa can be 
set 0, 0, 0] ■ This makes the pp system particularly effec- 
tive for measuring a. 

In B ^ pp decays the final state is a superposi- 
tion of CP-odd and CP-even states. An isospin-triangle 
relation Q holds for each of the three helicity ampli- 
tudes 74;^=_l o,+i, which can be separated through an 
angular analysis. The longitudinal polarization fraction 
/l = l^oP/ X] I^Ap can be determined through the mea- 
surement of the distribution of the helicity angles Oi and 
02, defined as the angles between the direction of the 7r+ 
and the direction of the B meson in the rest system of 
each of the tt^tt^ candidates. 

In this paper, we update our previous measurement Q 
of the branching fraction B and longitudinal polariza- 
tion fraction fi, in B^ — > p^p^ decays, along with B for 
P" — > p°/o, /o/oj /0°7r+7r~, and 7r+7r~7r+7r~. In addi- 
tion, we present the first study of the time-dependent 
CP asymmetry Acp in this mode. We determine the co- 
efficients C£° and S'£° of Acp for the longitudinal com- 
ponent, expressed as a function of At, the proper time 
difference between the two B decays in T(4S') — + B'^B^: 

Acp{At) = ~Cf cos AmAt -I- sin AmAt . (1) 

where Am = (0.507± 0.005)?i ps~^ is the mass difference 
between two P° mass eigenstates When combined 
with measurements of P" ^ r,+ n- cnri r+ 
Acp, B, and /l in P^ ^ 
analysis and improve the constraints on the penguin con- 
tribution to B pp decays. Changes with respect to 
our previous analysis [3| include a larger data sample, 
improved track-selection techniques, and inclusion of the 
P-decay time information. 

We use a sample of (465 ± 5) x 10^ 7(45") deca^ 
into PP pairs collected with the PaPAR detector 11 1 
at the PEP-II asymmetric-energy e+e~ collider [l3|. A 



and P+ ^ p+pO, 



p'^p'^ allow a complete isospin 



detailed description of the PaPAR detector is available 
elsewhere [1, [ll|. 

We select P M1M2 (7r+7r^)(7r+7r~) candidates, 
where Mi_2 stands for a p^ or /o(980) candidate, from 
neutral combinations of four charged tracks that are con- 
sistent with originating from a single vertex near the 
e^e" interaction point. We veto tracks that are pos- 
itively identified as kaons or electrons. The identifica- 
tion of signal P candidates is based on several kinematic 
variables. The beam-energy-substituted mass toes = 



[(./2-Hp.-PB)Vi?f 



where the initial e+e four- 



momentum (Pi, Pi) and the P momentum pe are defined 
in the laboratory frame, is centered near the P mass with 
a resolution of 2.6 MeV/c^ for signal candidates. The dif- 
ference AP = P™ — 1/3/2 between the reconstructed P 
energy in the center of mass (cm.) frame and its known 
value \/s/'2. has a maximum near zero with a resolution of 
20 MeV for signal events. Four other kinematic variables 
describe two possible tt+tt^ pairs: invariant masses mi, 
m2 and helicity angles 0i, 62- 

We use the kinematic selection of signal candidates de- 
scribed in We require 5.245 < mEs < 5.290 GeV/c^, 
|AP| < 85 MeV, 0.55 < mi,2 < 1-050 GeV/c^, and 
|cos0i 2| < 0.98. The extended di-pion invariant mass 
range allows us to study the non-resonant decays P'^ — *■ 
p'^TT+TT^ and B^ 7r+7r~7r+7r~, as well as P" — > p^ Jq 
and B'^ — > /o/o- The contributions from the higher mass 
resonances in this range are relatively small. We sup- 



press the dominant 



qq (q = u, d, s, c) continuum 



background using a neural network-based discriminant £, 
which combines eight topological variables |4] . 

We use multivariate P-flavor tagging algorithms 
trained to identify primary leptons, kaons, soft pions, 
and high-momentum charged particles from the other 
P, called Ptag [lH- The effective tagging efficiency is 
(31.1 ± 0.3)%. Additional background discrimination 
power arises from the difference between the tagging effi- 
ciencies for signal and background in seven tagging cate- 
gories (ctag). We determine At and its error a At from the 
spatial separation between the decay vertices of the signal 
P and Ptag and require \At\ < 15 ps and cta* < 2.5 ps. 

After application of all selection criteria, 72154 events 
are retained. On average, each selected event has 1.05 sig- 
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p°p° 



nal candidates, while in Monte Carlo (MC) [3, 
pies of longitudinally (transversely) polarized B'^ 
decays we find 1.15 (1.03) candidates. When more than 
one candidate is present in the same event, the candi- 
date that yields the smallest for the four-pion vertex 
is selected. Simulation shows that 18% of longitudinally 
and 4% of transversely polarized P" p'^ events are 
misreconstructed with one or more tracks that do not 
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originate from the B° p^p^ decay. These are mostly 
low-momentum tracks from the other B meson in the 
event. Such partially reconstructed candidates are in- 
cluded in the definition of the signal probability density 
functions (PDFs). 

We use an unbinned extended maximum likelihood fit 
to extract the B° p°p° event yield, /l, Cf , and 5*^ . 
We also fit for the event yields of B^ P^fo, B'~' 
/o/o, p^n+TT-, and B" tt+tt-tt+t:- decays. The 

likelihood function includes the background components 
from non-signal B decays and continuum. The PDFs for 
each component depend on ten discriminating variables: 
TTiEs, Ai?, f ,mi,m2,cos0i,cos02,ctag, At, and cta*- 

Since the statistical correlations among the variables 
are found to be negligibly small, we take the total PDF 
as the product of the PDFs for the separate variables. 
Exceptions are the kinematic correlation between the two 
helicity angles in signal, and mass-helicity correlations in 
several i?-decay classes and misreconstructed signal. 

We use double-Gaussian functions to parameterize the 
mES and AE PDFs for signal, and relativistic Breit- 
Wigner functions for the resonance line-shapes of and 
/o(980), with the /o(980) mass and width taken from [l^. 
The angular distribution for B decays [3] (expressed 
as a function of for B^ P^) is multiplied by 

a detector acceptance function determined from simu- 
lations. We assume that the p^ in B^ p°7r+7r~ is 
longitudinally-polarized (i.e., tt+tt" are produced in the 
S-wave), and we use the phase-space distributions for 
TT+vr^Tr+vr^. The (tttt) invariant mass and he- 
licity distributions of misreconstructed signal events are 
parameterized with empirical shapes in a way similar to 
that used for B background discussed below. The neural 
network discriminant £ is described by the sum of three 
asymmetric Gaussian functions with different parameters 
for signal and background distributions. 

The PDFs for non-signal i?-decay modes are generally 
modeled with empirical analytical distributions. Several 
variables have distributions identical to those for signal, 
such as Toes when all four tracks come from the same 
B, or TOi 2 when both tracks come from a p'^ meson. 
Also for some of the modes the two tt+tt^ pairs can have 
different mass and helicity distributions, e.g., when only 
one of the two combinations comes from a genuine p° or 
/o meson, or when one of the two pairs contains a high- 
momentum pion (as in B oitt). In such cases, we use a 
four- variable correlated mass-helicity PDF. The proper- 
time distribution for signal and background B decays is 
convolved with a resolution function while the time 
distribution of continuum background is assumed to have 
zero lifetime. 

The signal and S-background PDF parameters are 
extracted from simulation. The MC parameters for 
the Toes I AE, and £ PDFs are adjusted by compar- 
ing data and simulation in control channels with simi- 
lar kinematics and topology, such as -B" D~tt^ with 
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FIG. 1: Projections of the multidimensional fit onto the (a) 
niES, and (b) di-pion invariant mass m,r7r (average of mi and 
m2 distributions is shown), after a requirement on the signal- 
to-background probabifity ratio with the plotted variable ex- 
cluded, which enhances the fraction of signal events in the 
sample. This selection has 39% (60%) efficiency for signal for 
the rriES (m,r7r) projection. The data points are overlaid by 
the fuU PDF projection (solid black line). Also shown are the 
B" p'^p'^ PDF component (dotted fine) and the sum of all 
other PDFs (dashed fine). 

A'+TT^TT^. The continuum background PDF 
shapes are extracted from on-resonance sideband data 
(toes < 5.27GeV/c^), with parameters of most PDFs 
(for TOES, AE, £, 6*1, 02, and At) left free in the final fit. 
The tagging efficiencies, mistag fractions, and the param- 
eters of the proper-time distributions for signal modes are 
obtained in dedicated fits to events with identified exclu- 
sive B decays [13]. For inclusive B backgrounds these 
parameters are determined by MC and their systematic 
uncertainties are evaluated in data. 

We study the contributions of the dominant back- 
grounds by using high-statistics exclusive MC samples. 



We single out the i?'^ 



r=F, 5° -> p°K*° and B° 



JqK*^ modes, which have kinematic distributions simi- 
lar to those of the signal events, and parameterize their 
PDFs individually. The event yield for B^ a^n^ is al- 
lowed to vary in the fit, while the yields for B^ p'^K*'^, 
B^ foK*'^ are fixed to the expected values |17| . 

In addition, we construct a charmless event category 
consisting of B'^ P^P , B^ P^t^^, B^ P'^p'^, 
B+ a\TT+, B+ atfo, B rj'K, and B+ p°tt+ 
backgrounds. Kinematic distributions in these events, 
especially events in which at least one charged particle is 
not correctly associated to the B candidate, are similar to 
each other, and also to those of other, poorly measured, 
charmless decays. We allow the overall event yield for 
this category of events to vary in the fit, after fixing the 
relative weight of each mode to the expected value 
ITsj . The remaining events, which mostly originate from 
open 6 — > c transitions, are parameterized as a separate 
background component, with its yield left free in the fit. 

Table U summarizes the results of the fit. The B'^ 
p^ p^ decay is observed with a significance of 3.1 stan- 
dard deviations {a), as determined by the quantity S — 
■\/— 2 ln(£o/'Cmax), where £max is the maximum likeli- 
hood value, and is the likelihood for a fit with the sig- 
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TABLE I: Event yields; Iraction of longitudinal polarization (/l); selection efficiency (Eff ) corresponding to measured polariza- 
tion; branching fractions (B); branching fraction upper limits (UL) at 90% confidence level (CL); and significance S including 
systematic uncertainties. First uncertainty is statistical and second is systematic. 



Mode 
















fa fa 

















Yield /l Eff (%) B (IQ-'') UL (10~^) S (a) 

99l^^± 15 0.75t;5 i4 ± 0.04 23.28 ± 0.07 0.92 ± 0.32 ± 0.14 - 3.1 
3l2o±5 - 24.16 ±0.09 0.03j:E;:i^ ± 0.05 < 0.34 
+**± 2 - 27.22 ±0.07 0.05+^'^'? ± 0.02 < 0.16 



1.68 ±0.01 -i.allj 

0.55 ±0.01 3.2t^'^5'^ ± 3.4 < 21.1 



-12t^^±9 - 1.68 ±0.01 -1.21^-9 ±1.1 < 8.7 



-35 
5+30 



nal contribution set to zero. Both likelihoods include sys- 
tematic uncertainties, which are assumed to be Gaussian- 
distributed and are discussed below. This significance 
level corresponds to a probability of 1.0 x 10"'^ that the 
observed signal yield is consistent with a background 
fluctuation. We do not observe significant event yields 
for 5° p°fo or B° /o/o decays, nor for the non- 
resonant decays p°7r+7r~ and i?" 7r+7r^7r+7r^. 
We find 280 ± 53, 670 ± 96, 2329 ± 147, and 68701 ± 281 
events for the — > afir^, charmless, b ^ c, and contin- 
uum backgrounds, respectively, consistent with expecta- 
tions. In Fig. [1] we show the projections of the fit results 
onto mEs and m^jr- 

From the fit to the proper-time distribution of the data 
sample, we determine the CP-violating parameters 



0.3 ±0.7 (stat.)±0.2 (syst.) 
0.2 ±0.8 (stat.)±0.3 (syst.) 



for the longitudinal component of the — > p'^p^ sam- 
ple. The statistical correlations between fit parameters 
for i?° p'^p'^ are given in Table [TTl 



TABLE II: Correlation matrix for B° 



(P (P parameters. 



Parameter 


Yield 


fL 




cf 


Yield 


1.000 


0.086 


-0.136 


-0.273 


fL 




1.000 


-0.006 


-0.174 


sf 






1.000 


-0.035 


cf 








1.000 



Dominant systematic uncertainties in the fit originate 
from statistical errors in the PDF parameterizations due 
to the limited number of events in the control samples, 
variations in the B background branching ratios fixed in 
the fit, and from potential fit bias. The PDF param- 
eters are varied by their respective uncertainties to de- 
rive the corresponding systematic errors. The fit bias is 
studied in a large number 1000) of MC experiments, 
in which signal and dominant charmless B background 
events are fully simulated, while other backgrounds are 
sampled from their respective PDFs. We correct for the 
bias of 7.9 ± 2.0 events for p^ p^ , 0.03 ± 0.02 for 



and 0.07 ± 0.03 for 5*^°. The uncertainties associ- 
ated with the B background model are 4 events for the 
signal yield, 0.01 for /l, 0.01 for Cf and 0.11 for Sf. 
The systematic uncertainties due to the charmless back- 
ground composition, arising from the uncertainties in the 
individual branching ratios and the CP content of the B 
background 18, 3|, are 5 events for the signal yield, 0.01 
for /l, 0.18 for and 0.14 for Sf. The above system- 
atic uncertainties do not scale with event yield and are 
included in the calculation of the significance of the re- 
sult. 

We estimate the systematic uncertainty due to the in- 
terference between the B'^ p^ p^ and aj'^7r± decays us- 
ing simulated samples in which the decay amplitudes for 
B^ — > p^ p^ are generated according to this measurement 
and those for 5° o^tt^ correspond to a branching 
fraction of (33.2 ± 4.8) x lO"*^ The strong phases 

and CP content of the interfering o^'k^ state are varied 
between zero and a maximum value using uniform distri- 
butions. We take the RMS variation of the fitted values 
(13 events for the p"p" yield, 0.02 for /i, and 0.04 for 
S"™ and C™) as a systematic uncertainty. 

Uncertainties in the reconstruction efficiency arise from 
track finding and particle identification, and are deter- 
mined by dedicated studies on control data samples. Un- 
certainties due to other selection requirements, such as 
vertex probability, track multiplicity, and thrust angle, 
amount to 2.4% for the event yields, and are negligible 
for the polarization and CP observables. 

We perform an isospin analysis oi B ^ pp decays, 
by minimizing a that includes the measured quanti- 
ties expressed as the lengths of the sides of the isospin 
triangles We use the measured branching fraction 
and fraction of longitudinal polarization of S"*" 
decays Q, the measured branching fraction, polariza- 
tion, and CP parameters S'^~ and determined 
from the time evolution of the longitudinally polarized 
B^ p^p^ decays 0, and the branching fraction, po- 
larization, and CP parameters Sj^ and of B'^ p^ p^ 
reported here. We assume uncertainties to be Gaus- 
sian distributed and neglect 7 = 1 isospin contributions, 
electroweak loop amplitudes, non-resonant, and isospin- 
breaking effects. 
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FIG. 2: Confidence level (CL) on a — a^g determined from 
the isospin analysis. The long-dashed curve is obtained with- 
out the two CP parameters S'l^ and C£''. The dotted curve 
corresponds to the isospin analysis without Sj^ , and the solid 
curve CL includes both C£'' and in the fit. The horizontal 
dotted lines correspond to the 68% (top) and 90% (bottom) 
CL intervals. 



We obtain a 
15.7° (< 17.6°) 



68% (90%) CL limit \a - a^sl < 
where a^ff is defined by sin(2acfi) = 
Fig. [2] shows the confidence level 
in the isospin anal- 
ysis fit. We observe four solutions near zero, as in the 
_B — > TTTT isospin analysis 2l|. The additional constraint 



with and without use of Sj^ and 



from 



provides some discrimination among the four 
solutions. 

In summary, we confirm our earlier evidence for 
j^Q _^ pOpO (jgcays with a significance of 3.1a and mea- 
sure the branching fraction, longitudinal polarization 
fraction, and CP asymmetries in these decays. These 
measurements combined with those for _B+ P'^ and 
p^p~ decays provide a constraint on the pen- 
guin uncertainty in the determination of the CKM uni- 
tarity angle a. We find no significant evidence for the 
decays B° ^ pVo, B° ^ /o/o, B° ^ pO^+^- or 
B^ — > Tr+vr^Tr+vr^. 
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